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The Schrodinger motion of a charged quantum particle in an electromagnetic potential can be 
simulated by the paraxial dynamics of photons propagating through a spatially inhomogeneous 
medium. The inhomogeneity induces geometric effects that generate an artificial vector potential to 
which signal photons are coupled. This phenomenon can be implemented with slow light propagating 
through an a gas of double-A atoms in an electromagnetically-induced transparency setting with 
spatially varied control fields. It can lead to a reduced dispersion of signal photons and a topological 
phase shift of Aharonov-Bohm type. 



Introduction. — It is known since the ground-breaking 
work of Berry on geometric phases [H that artificial gauge 
potentials can be induced if the spatial dynamics of a sys- 
tem that obeys a wave equation is confined in a certain 
way. For instance, if the internal Hamiltonian of neu- 
tral atoms contains an energy barrier but the spin eigen- 
states are spatially varying, gauge field dynamics can be 
induced 0]. In the limit of ray optics, moving atomic en- 
sembles could simulate the propagation of light around 
a black hole or generate topological phase factors of the 
Aharonov-Bohm type 0], and inhomogeneous dielectric 
media could generally exhibit geometric effects such as 
an optical spin-Hall effect and the optical Magnus force 

a. 

In this paper, we propose to use electromagnctically in- 
duced transparency (EIT) to generate an artificial vector 
potential for the paraxial dynamics of signal photons that 
simulates quantum dynamics of charged particles in a 
static electromagnetic field. Not only the ray of light but 
also its mode structure is affected, resulting in a paraxial 
wave equation that is equivalent to the Schrodinger equa- 
tion for charged particles. Furthermore, the form of the 
artificial vector potential can be easily controlled through 
spatial variations in the control fields. We suggest con- 
figurations that generate homogeneous quasi-electric and 
magnetic fields as well as a vector potential of Aharonov- 
Bohm type. 

Although the treatment in this paper is based on EIT, 
the effect presented here is more general: it will occur in 
any medium that supports a set of discrete eigenmodes 
for a propagating signal fields with different indices of 
refraction. If the parameters governing these eigenmodes 
vary in space, the signal modes will adiabatically follow, 
acquiring geometric phases that affect their paraxial dy- 
namics. 

Review of EIT with multi-A atoms. — The effect takes 
place in an atomic multi-A system, in which two ground 
states are coupled to Q excited states by Q pairs of con- 
trol {Vt q ) and signal (a q ) fields (Fig. [T]). An experimen- 
tally relevant example of such system is the fundamental 
Dl transition in atomic rubidium, where both the ground 
and excited levels are split into two hyperfine sublevels 
[H . We assume that the detunings are small so each sig- 



nal field a q interacts only with the respective transition 
\B) <-» \A q ) with the associated atomic operator (?B,A q 
and vacuum Rabi frequency g q . In this case, the parax- 
ial wave equation for each signal mode can be cast into 
the form 

/ ic \ 

where the wave propagates along the z axis, Aj_ = d% + 
dy, N is the number of atoms and k is the wavevector 
which we assume approximately independent of q. In 
Ref. Q we have constructed a unitary transformation 

Q 

a q = J2 W v Is (2) 

s=l 

that maps the original field modes a q to a new set of 
modes b qi such that one and only one of the new modes, 

Q 

9=1 

(where R q = il q /(g q il±) and fij_ = \jYl%i \^ql9q\ 2 de- 
pend on the control fields) couples only to an atomic dark 
state and experiences EIT [1, 0, 0] • All other superposi- 
tions of field modes are absorbed. This transformation is 
given explicitly by W qq > = juiqW*, — 5 qq > , with 7 = Rq + 1 
and w q = r y~ 1 (S Qq + R q ). 

The EIT mode 6q interacts with the multi-A atoms 
in the same fashion as does the signal field in a regu- 
lar 3-level system. While propagating through the EIT 
medium, it gives rise to a dark-state polariton associ- 
ated with zero interaction energy Q. All other modes 
couple to atomic states whose energy levels are Stark- 
shifted by the interaction with either the pump field or 
the other signal modes b q (q 7^ Q). The resulting energy 
gap guarantees that, if the amplitudes and phases of the 
control fields are slowly changed, the composition of the 
dark-state polariton, and hence the EIT mode 6q, will 
adiabatically follow. It has been proposed and ex- 
perimentally demonstrated [H that a variation in time of 
the control fields can therefore be used to adiabatically 
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FIG. 1: Multi A-system: Q excited states \A q ) are each cou- 
pled by a classical control field Q q to the ground state |C) 
and by a quantized field a q with detuning 8 to state \B). 



transfer optical states between signal modes. In this pa- 
per, we focus on spatial propagation of the EIT mode 
under control fields that are constant in time, but varied 
in space. 

Derivation of the gauge potential. — We proceed by ex- 
pressing Eq. |T]) in terms of the new signal modes b q . 
Employing the vector notation a = {ai,--- ,clq} and 
ob,a = {giVB.Ai > • • • , 9QVB,A Q } we get 

di + 4 z -Tk A -W b = iN ^- < 4 > 

Throughout the paper, the double arrow denotes aQx 
Q matrix. Because W depends on space and time, the 
differential operators have to be applied to both W and 

b. As a result, transformation @ brings about additional 
terms into the equation of motion, that can be written 
in form of a minimal coupling scheme by introducing the 
Hermitian gauge field 



A = iW^diW, 



(5) 



where i = t,x, y, z. We multiply both sides of Eq. (|U) by 
W> and exploit the unitarity of W to show that diW' = 
~W\diW)W ] from which it follows that -W^dfW = 
A\ + idiAi. The dynamic equation for the b modes can 
then be written as 



id t + A^ b = - (icd z + cA z ^ b 



(6) 



+77rHv ± - A^fb - w^Nvba 

Ik 

with Vj_ = (d x ,d y ). This equation has the structure of 
a 2+2 dimensional field theory with minimal coupling. 

Under the assumption that the control fields do not 
depend on t and z we can make a temporal Fourier trans- 
formation of the slowly varying amplitudes, which results 
in the paraxial wave equation 



id z b(S) 



- 6 - + ±:(-iv±-A ± y 

^ + A 
-W^-a BA {5). 

c 



The gauge potential is given explicitly by 

«-. Q 

A± = i R* ( Vj_ R q ) vSvft — ry(V \w)vft +i"f*wV ^w^ . 

q=l 

(8) 

The full matrix Aj_ is a pure gauge: it has emerged solely 
as a consequence of the unitary transformation ([2]), which 
reflects our choice to describe the system in terms of the 
new modes b q rather than the original modes a q . How- 
ever, this choice is motivated by the fact that the EIT 
mode 6q is the only mode that is not absorbed. Ab- 
sorption of other modes b q (with q ^ Q) means that 
the index of refraction for these modes has a significant 
imaginary part. This separates the EIT mode 6q from 
other 6-modes and ensures that it will adiabatically follow 
variations of the control fields. Therefore, when analyz- 
ing the evolution of 6q, we can neglect the off-diagonal 

terms in the matrix (— iVj_ — A±) 2 in Eq. ([7]) and write 



id z b Q (5) 



(9) 



-(W^ff BA ) Q {6)- 6 -b Q (6) 
-7^HV ± - A±) Qq (-iV± - A ± ) qQ b Q {S). 



This equation does not include the whole matrix Aj_. 
Consequently, this potential no longer acts like a pure 
gauge but attains physical significance in determining the 
spatial dynamics of the EIT mode. 

The first term on the right-hand side of Eq. @, re- 
sponsible for the interaction of the light field with the 
EIT medium, takes the same form as the susceptibility 
of EIT in a single A-systcm. Neglecting decoherence, we 
can write it as @ (W^ctba)q(o~) = ^~^Q: wltn tnc 
EIT group velocity veit = c£l 2 /{Ng 2 ). Note that weit 
depends on the spatial position because Q, docs. This 
transforms Eq. §§§ to 

i a.S =a ( - i Vx-W-^ + |)8« do) 



with 



Aqq = iJ2R* q V±R q = -J2\R q \ 2 V±Arg(R q ), 



q=l 
3-1 



9=1 



(7) 



$ = £K^)q/ = -^w + EI v ^I 2 (ii) 

9 =1 q=l 

being, respectively, the "quasi-vector" and "quasi-scalar" 
potentials. 

We see that the paraxial spatial evolution of the EIT 
signal mode is governed by the equation that is identi- 
cal (up to coefficients) to the Schrodinger equation of a 
charged particle in an electromagnetic field. This is the 
main result of this work. By arranging the control field in 
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a certain configuration, one can control the spatial prop- 
agation of the signal mode through the EIT medium. 

Some steering of the EIT mode is possible even in a 
single- A system by affecting the term S/veit in Eq. (fl"0|) . 
which results in nonuniform refraction for this mode 
ficl The action of quasi-gauge fields (fTTjl is fun- 

damentally different: deflection of the signal field occurs 
not due to refraction (the refraction index on resonance 
is 1), but due to adiabatic following. 

The case of two control fields: homogeneous electric 
and magnetic quasi-fields. Of particular practical im- 
portance is the simplest non-trivial case with Q = 2. 
We parametrize the control fields by writing Rx,2 — 
y/l/2 ± Re l<y ^ ±e \ The corresponding Rabi frequencies 
are then £li = h(x,y) giRi, with h(x,y) being an arbi- 
trary common prcfactor. This paramctrization yields the 
gauge potentials 



= -V_l0 - 2R\7 ± 6 
(V ± i?) 2 



(12) 



1-4R 2 



(V±0) 2 (l-4i? 2 ). 



Similarly to usual electrodynamics we can use a gauge 
transformation [lij], Aqq — Aqq + Vj_/, to eliminate 
the term Vj_0 from Eq. (|T2]) . The common phase (f> of 
the control fields therefore does not contribute and can 
be set to zero. 

A simple way to generate a term that corresponds to a 
one-dimensional scalar potential V(x) for a Schrodingcr 
particle is to choose R = and = f Q dx' y / 2kV(x') 
This choice of control fields leads to Aqq = and $ = 
2kV{x). 

For the special case of a constant electric quasi-field 
along the x axis, V(x) = —Fx and subsequently 



6 = -^4fcF|x| 3 /3, 



(13) 



where x < is assumed for the region of interest. A res- 
onant (S = 0) Gaussian solution to Eq. (fTO]) is displayed 
in Fig. [21a). The center of the Gaussian beam is shifted 
by an amount x c tr = Fz 2 /2k, which is equivalent to the 
motion of a charged particle in a constant electric field. 

The control field phase profile (fT3"|) can be implemented 
using, for example, a phase plate. The assumption that 
the control fields do not depend on z implies that the 
Fresnel number for these fields must be above 1, i.e. that 
the characteristic transverse distance over which these 
fields significantly change must be larger than ~ V XL, 
where L is the EIT cell length. This imposes a limi- 
tation on the magnitude of the electric quasi-field: from 
Eq. CE3]) we find F < A" 1 / 2 ^ 3 / 2 and thus x ctI < sfhjX. 
Assuming that the signal field also has a Fresnel number 
of at least 1, and thus satisfies 2zr > L (with Rayleigh 
length zr = kw 2 /2, w being the signal beam width at 
the cell entrance), we find that in a realistic experiment, 
the maximum possible signal beam displacement due to 



the quasi-electric field is on the order of the signal beam 
width w. 

To generate a homogeneous magnetic quasi-field along 
the z-axis the quantity B = V x Aqq = 2V ±9 x V^i? 
should be constant. However, it seems difficult to si- 
multaneously achieve a vanishing electric quasi-field E = 
— Vj_$. A choice that minimizes the electric quasi- 
field around the origin is given by 9 = ^~B]2x and 
R = \J Bj2 y. The quasi-potentials then become Aqq = 
—B y e x , which corresponds to the Landau gauge in stan- 
dard electrodynamics, and <f> = B + 2B^y A + 0(y 6 ). If 
$ is neglected, a Gaussian solution to the paraxial wave 
equation is given by 

oq = Afcscu(z) exp cotu(z)Ax 2 + Ax • p c 

-§AxAy - jgPcxPc.y + \VcPc,y)\ (14) 

where we have set Ax = (x — x c ,y — y c ), u(z) = 
Bz/(2k) - itanh _1 (277) and n = Bw 2 /A. Here x c = 
xo + {k/B)(x' s'm(Bz/k) + Xq(1 — cos(Bz/k)) denotes 
the classical spiral trajectory of a charged particle in a 
magnetic field, with x^, = dx c /dz, initial position xo and 
initial velocity Xq. For convenience we also have defined 
Xq = (j/q, — a^o) an d the classical canonical momentum 
p c . We remark that p c ^ x is a constant of motion. The 
evolution of the signal mode is displayed in Fig. 2(b). 

A surprising feature of solution (|T4|) is that the diffrac- 
tive divergence of the signal beam is reduced: the width 
squared of the Gaussian, 

2 w 2 ( , , „ ,Bz s 

— ; 7 = 1 + 4?/ 2 - (1 - 4n 2 ) cos — ) 

Re(iBcotK) V V k 

(15) 

varies periodically with z instead of monotonically in- 
creasing. This effect is known for electron wavepackets 
fl3| and can be understood as a consequence of the cir- 
cular motion of particles in a magnetic field: instead of 
dispersing, two-dimensional particles in a magnetic field 
will simply move on circles of different size (depending on 
their velocity), but with the same angular velocity. The 
particle cloud will therefore not spread but "breathe" . 

It remains to show that non-adiabatic coupling to other 
modes can be suppressed for realistic experimental pa- 
rameters. This is the case if the strength of the gauge 
field terms coupling &q to other modes, which for the 
quasi-magnetic field are of the order B/(2k), are much 
smaller than the difference in the respective linear sus- 
ceptibilities xi- F° r the EIT mode 6q, xi = <V^eit 
with veit defined above Eq. (flQ|) : for the other modes 
it can be approximated by the susceptibility of a two- 
level medium, xi = —^Ng 2 {5 — i^/2)/(c^ 2 ). Evalu- 
ating this relation at resonance leads to the condition 
?/ <C (fcw) 2 n37r/2, with n = N/(Vk 3 ) being the number 
of atoms in the volume k~ 3 , which can easily be fulfilled 
in an experiment. 

Aharonov-Bohm potential for photons. — One of the 
most intriguing phenomena of charged quantum particles 
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FIG. 2: Paraxial propagation of a signal beam over twice the 
Rayleigh length in the presence (solid) and absence (grey) of 
a constant (a) electric field along the x axis and (b) magnetic 
field along z. The dashed line represents the center of the 
grey beam. The effect of the fields is somewhat exaggerated. 

in electromagnetic fields is the Aharonov-Bohm (AB) ef- 
fect [IH • Its two astonishing features are (i) a phase shift 
induced by the vector potential in a region in which elec- 
tric and magnetic fields are absent, and (ii) its topological 
nature: the phase shift does not depend on the particle 
trajectory as long as it encloses a magnetic flux. Because 
(unlike genuine electromagnetism) the potential (|5|) is a 
differential function of the control fields, it is impossible 
to simulate feature (i) with quasi-charged photons. How- 
ever, we will show here that a mathematically equivalent 
topological phase shift does exist for the optical case. 

To generate an AB potential for photons we propose 
to use two counter-rotating Laguerre-Gaussian control 
fields, i.e., fields that possess an orbital angular momen- 
tum. If these control fields are spatially wider than the 
signal fields, the corresponding Rabi frequencies can be 
approximated in cylindrical coordinates (r, ip) by Qi = 
giSire lLp and Q 2 = gisire~ %v . The gauge potentials (fl~2"l) 
then become A QQ = -2R/r e v and $ = (1 - 4i? 2 )/r 2 , 
with R = i(| Sl | 2 - |s 2 | 2 )/(|si| 2 + |s 2 | 2 ). The potential 
Aqq corresponds exactly to an Aharonov-Bohm poten- 
tial for charged particles as it is created by a solenoid. 

Solutions of the paraxial wave equation (fTU)) can be 
found in cylindrical coordinates by expanding the field 
mode as 6q = r -1 / 2 dime's, B m (z, r) exp(imtp). Because 
of Q ~ r, the EIT group velocity can be written as ^eit = 
vr 2 with v = c-\/|si| 2 + |s 2 | 2 /A. Exact solutions are 
given by Bessel functions, 

B rn = e-- 2 ^ 2 *) (a m ^J v {nr) + p m \f^Y v { K r)) , 

(16) 

with v = vT-j- m 2 + ARm — 2kvS. For monochromatic 
signal fields this corresponds to a rotation of the trans- 
verse mode structure. For R = ±1/2 the potential trans- 
fers a unit amount of angular momentum to the signal 
light, but generally the amount can vary continuously be- 
tween —h and h. Signal photons in the EIT mode there- 
fore form a two-dimensional bosonic quantum system in 
an Aharonov-Bohm potential. 

Conclusion. — We showed that EIT in a multi-A sys- 
tem can be used to generate a variety of geometric ef- 



fects on propagating signal pulses that mimic the be- 
havior of a charged particle in an electromagnetic field. 
We found specific arrangements of two spatially inhomo- 
geneous pump fields in a doublc-A system which gener- 
ate quasi-gauge potentials which correspond to constant 
electric and magnetic fields. Furthermore topological ef- 
fects like the Aharonov-Bohm phase shift can be induced. 
The latter is significantly different from the proposal of 
Ref. Q in that it is based on spatially inhomogeneous 
pump fields rather than the Doppler effect in moving me- 
dia. 

This paper investigated EIT in systems with two 
ground levels. In such a system, there is only one EIT 
mode, which results in an Abelian U(l) gauge theory, 
making the physics analogous to electromagnetism. By 
extending to multiple ground levels, it may be possible to 
obtain multiple EIT modes and model non- Abelian gauge 
potentials. This will be explored in a future publication. 
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